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The muon transfer probabilities between muonic hydrogen and an oxygen atom are calculated 
in a constrained geometry one dimensional model for collision energies between 10~® and 10^ eV. 
For relative translational energies below 10~^ eV, for which the de Broglie wavelength (> 1 A) is 
much larger than the characteristic distance of the potential interaction (~ 0.1 A), the problem 
corresponds to an ultra-cold collision. The close-coupling time-independent quantum equations 
are written in terms of hyperspherical coordinates and a diabatic- by-sectors basis set. The muon 
transfer probabilities are qualitatively interpreted in terms of a model involving two Landau-Zener 
crossings together with the threshold energy dependence. Based on this analysis a simple procedure 
to estimate the energy dependence of the muon transfer rate in three dimensions, is proposed. These 
estimated rates are discussed in the light of previous model calculations and available experimental 
data for this process. It is concluded that the high transfer rates at epithermal energies inferred 
from experiments are unlikely to be correct. 



I. INTRODUCTION 

Negative muon transfer between exotic atoms (muonic 
hydrogen, for instance) and other atoms or molecules, 
has been extensively studied in the framework of muon 
catalysed nuclear fusion (see Ref. and literature cited 
therein). Also, the structural and spectroscopic proper- 
ties of these species is of interest for metrology as well as 
a test of quantum electrodynamic theories j3| . 

Very recently, a new method to measure the hyperfine 
structure {F = 0, 1) of muonic hydrogen (p^) based on 
the coUisional energy dependence of muon transfer from 
the muonic hydrogen to an oxygen molecule, has been 
proposed When muonic hydrogen in the -F = hy- 
perfine ground state is laser excited to the F = 1 level, 
collisions with H2 convert this excess energy into kinetic 
energy, giving an additional 0.12 eV translational energy 
to the muonic hydrogen. When the muon is transfered 
to an oxygen atom, it is captured in high (n — 5, 6) 
states which promptly de-excite and emit X-rays. If the 
muon transfer rate to oxygen changes significantly from 
thermal (0.04 eV) to epithermal (0.16 eV) energies, the 
measurement of the X-ray emission intensity with and 
without laser excitation can be used for the determina- 
tion of the hyperfine splitting in muonic hydrogen 0, Q . 
This proposal was based on the work of Werthmiiller et 
at Q which indicates that muon transfer from hydrogen 
to oxygen increases by a factor of 3-4 going from thermal 
to epithermal collision energies. 

While the experimental transfer rate at thermal ener- 
gies have been borne out by the calculations of Sultanov 
and Adhikari 0, the rates at epithermal energies have 
never been confirmed theoretically. Even more, such an 
increase of the transfer rate with energy does not have 
a clear theoretical explanation. Werthmiiller et al |J| 
suggest the existence of a resonance at low energies. It 



is worth to stress that such an increase has not been 
observed in the case of muonic hydrogen colliding with 
sulfur while for the case of CH4, experimental data 
suggest that the transfer rate actually decreases going 
from thermal to epithermal energies Q- In order to im- 
prove our understanding of thisprocess and to assess the 
validity of the proposal of Ref. 3] , it is important to per- 
form model calculations to study in detail this reaction 
as a function of energy. It is the purpose of this paper to 
present such a study. 

Since the muonic hydrogen has to approach one of the 
oxygen nuclei very close in order for the muon to transfer 
S|, the process can be described as 
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Although there have been several full three- 
dimensional calculations of muon transfer rates at 
low energies between muonic-hydrogen and low-Z atoms 
(see literature cited in Ref. Q), there is none when 
the transfer involves nuclei with Z > Z . There are 
several good reasons for that. As Z increases (Z = 8 in 
our case) there is a larger initial-channel polarization 
and a stronger final-channels Coulomb interaction. In 
addition the number of open channels even at zero 
relative kinetic energy increases with Z making the full 
three-dimensional calculation very heavy. Thus up-to- 
now only approximate calculations have been performed 
for the muon transfer rate between muonic-hydrogen 
and oxygen |^ Among them, the most recent 

and accurate is the one done by Sultanov and Adhikari 
who used a two-channel approximation to the 
integro-differential Fadeev-Hahn formalism to calculate 
the muon transfer rates from muonic-hydrogen to the 
n = 5 states of oxygen. 

We have chosen to perform exact calculations but in 
a restricted configuration with the muon moving along 
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the line joining the proton to the oxygen atom. Indeed, 
for this system the cohnear configuration is the most fa- 
vorable for the transfer process. Furthermore, since the 
colinear configuration provides the minimum energy path 
for the reactance channel, even for an initial non colinear 
configuration there will be efficient orientational effects 
in particular for low translational energies. 

Two model interaction potentials have been used in 
the calculations: 

1. A pure coulombic potential with the bare oxygen 
nucleus O*^. 

2. A shorter range model potential with a distance 
dependent effective charge on the oxygen nucleus 
(calculated via a Thomas-Fermi electronic density 
model) to simulate the effect of the outer electrons. 

The Hamiltonian was written in terms of hyperspher- 
ical coordinates. A piecewise diabatic basis set on the 
hyperspherical angle was used to expand the wave func- 
tion. The basis functions have in turn been expanded 
in terms of first order Legendre functions to efficiently 
handle the two coulombic singularities at fixed hyper- 
radii. The resulting close-coupling time-independent 
Schrodinger equations in the hyperradius were solved us- 
ing a de Vogelaere algorithm and the partial and total 
muon transfer probabilities were determined by the stan- 
dard S-matrix analysis at large distances. Since for en- 
ergies below 10 ~^ eV the muon transfer process studied 
here is equivalent to an ultra-cold collision (de Broglie 
wavelength (> 1 A) much larger than the characteristic 
distance 0.1 A) of the potential interaction), special 
care had to be taken to the asymptotic analysis in the 
reactant channel. 

The paper is organized as follows. Section 2 introduces 
the model and the methodology used in the calculations. 
Section 3 presents the calculated muon transfer probabil- 
ities together with their interpretation in terms of simple 
Landau-Zener and threshold models. A procedure to es- 
timate their energy behavior in the 3-dimensional case, 
is also presented. In section 4, the muon transfer rates 
are discussed in the light of previous model calculations 
and available experimental data for this process. Finally, 
section 5 is devoted to the conclusions. 
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FIG. 1: Thomas- Fermi effective charge for the oxygen atom 
as a function of the distance to the O*^ nuclei. 
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The hyperspherical coordinates are then defined by 



tan d = r/R 



(6) 



where < 6 < 9fj_ , with 0^ = arctan(TOp/m). In terms 
of these coordinates and after regularisation of the wave 
function by the factor ^/p (the volume element is then 
given by dpd9), the time-independent Schrodinger equa- 
tion at total energy E is: 
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where V{p, 9) is the interaction potential given by 



II. METHODOLOGY 



The mass-scaled Jacobi coordinates adapted to the en- 
trance channel of reaction |^ are : 
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Two sets of calculations have been performed. One 
(2) with the bare oxygen nucleus, for which = 8 in ((HJ, 
and another with (d) being an effective charge given by 
the Thomas-Fermi model ^| and represented in Fig. 
^ ' The total wave function il^ip, 9) is expanded in terms of 

basis set wave functions depending on the hyperspherical 
angle 9. We use a diabatic-by-sector representation. In 
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each sector p„ — Spn < p < Pn + Spn',n = 1, ...,iVp we 
write: 

^(p,0)=^i^,(p)(/.,(0;p„) (9) 

i 

where 4'i{9] p„) are eigenstates of the Hamiltonian at 
fixed Pn distances. Their calculation requires the solu- 
tion of a bound state problem for a potential showing 
Coulomb singularities aX 9 = and 6 = 9^. Treatment of 
Coulomb singularities requires the use of specific meth- 
ods such as increasing the density of grid points, or alter- 
natively the oscillation frequency of the basis functions, 
near the singularities. Mapped Fourier J^il^], Lagrange 
[T^ or Schwartz |0, interpolations are examples of 
such methods which have been used for cases with one 
singularity. 

In our case we have two singularities to deal with at 
fixed hyperradius. Therefore we have used the auxilliary 
coordinate x = 2 9/9^ — 1. This new variable is equivalent 
to one of the two hyperspherical elliptic angles [T^l, but 
constrained to be in the range [-1,1] with the singulari- 
ties being at the boundaries. If the </>i(0; Pn) functions are 
renormalized as: 4>i{9; p) = (1 — x'^)^ (f>i{x; p), both sin- 
gularities of the potential are regularized. The (j)i{x] p) 
functions obey the differential equation: 

( %^ D+il- x'){V -^)) P) 

\ m9f, p^ Smp-^ J 

= il-x^)e,ip)Mx;p) (10) 

where 

1,^(1-.^)^ -24-^ (11) 

Eq.^lis solved by expanding p) functions on the 
basis set of the eigenvectors of £>, which are the associ- 
ated Legendre functions P^{x). This yields the general- 
ized eigenvalue problem : 

\ m9f^p^ J 

= e,(p)O0,(p) (12) 

where 4>i{p) is the vector of the unknown coefficients of 
the function 4>i{x] p) in the Legendre basis set. The ma- 
trix representation of the differential operator D is di- 
agonal, with diagonal elements —n{n + 1). W is the 
potential coupling matrix and O is the representation of 
(1 — x'^) in the Legendre basis. These matrices are calcu- 
lated from a transformation of their diagonal representa- 
tion in a grid of Gauss-Legendre quadrature points [T8l |. 
The system ((T^ is transformed to a standard eigenvalue 
problem after left multiplication by 0~^. The method 
is expected to converge fast with the size of the basis set 



(or equivalently with the number of grid points) since the 
use of a Gauss-Legendre scheme provides an adequately 
high density of points near x = ±1. Since the wavefunc- 
tions concentrate in the vicinity of x — ±1 as energy 
decreases or as p increases, we expect convergence to be 
most difficult for lower levels and large p. At large p, we 
also expect convergence to be more difficult to achieve 
for Ofi states which are more compact than pp. states. 
We have computed the relative error for a given number 
riL of Legendre basis functions, using ul = 450 as a ref- 
erence. For riL — 150, we obtained 86 states converged 
better than 10~^ of relative error up to p ~ oq. For larger 
p values, the convergence starts to deteriorate, first for 
the ground state Op{n = 1), then for the first excited 
state n = 2 near p = 2 oq... Using larger values of de- 
lays the limit of convergence deterioration. For instance, 
with fiL = 250, the ground state is converged with 10"^ 
relative precision up to p ~ 5ao, up to p ~ 10 oq with 
riL = 350. This latter value of was used to generate 
the results shown below, and provided a basis converged 
better than 10~^ up to p ~ 40 oq for all states except 
the first four ones which do not play any significant role 
in the charge transfer process considered here. As ex- 
pected, among all these states, the ones corresponding 
to bound states of pp are the easiest to converge : the 
relative error at p = 40 gq for pp{n=l) is 10^^" with 
riL — 250 and 10^^^ with ul — 350. High accuracy 
on the pp{n = 1) channel is especially important for our 
present study since we consider initial relative kinetic en- 
ergies as low as 10~^ eV. 

The close-coupling equations are integrated along the 
hyperradius p using the de Vogelaere algorithm [T^ . 
In order to have convergence of the transfer probabil- 
ities to better than 1 % in the energy range consid- 
ered here, we included 29 channels: {{p p)n=i-4 + O and 
p + (/i 0)„=i-25)- The integration of the coupled equa- 
tions was performed from the origin to Pgjjfj ~ 30 oq. 

The asymptotic analysis has been performed using the 
appropriate Jacobi coordinates for the entrance and for 
the product channels. The asymptotic states were ex- 
pressed as products of Coulombic bound states and trans- 
lational functions. For the product channels, the trans- 
lational functions were approximated by simple plane 
waves. However, since the initial relative translational 
energy can be very small, the entrance channel has to 
be described with care to account for residual long range 
interactions. For the pure coulombic potential, this resid- 
ual interaction has a 1/i?^ dependence (charge-dipole in- 
teraction). In this case, Bessel functions of imaginary or- 
der have been used as a basis set in R. For the screened 
Thomas-Fermi potential, wave functions depending on 
R have been numerically computed by inward integra- 
tion from very large distances. These asymptotic wave 
functions and their derivatives are projected onto the hy- 
perspherical basis functions (j)i{x]pQ^^) [19^ . Transfer 
probabilities are then obtained by a standard S-matrix 
analysis. 
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FIG. 2: Adiabatic energies as a function of the hyperradius. 



1 



total 






n — 5 \ 








10"^ 10"^ 10"^ 0.01 0.1 


1 


10 100 100 



E (eV) 

FIG. 3: Partial transfer probabilities calculated with the pure 
coulombic potential (bare oxygen nucleus). 



III. RESULTS OF THE CALCULATIONS 

We have performed calculations for collision energies in 
the range 10~^— 10^ eV. In Fig.|21we present the adiabatic 
energies (see Ea. llO|) as a function of the hyperradius p. 
The origin of energies has been chosen to coincide with 
the asymptotic limit of the entrance channel (p/Lt)„=i+0. 
Thus the calculations cover the energy range between this 
limit and the p + (/i 0)„=io threshold. 

The dynamics of muon transfer can be qualitatively 
understood by inspection of this figure. Starting in chan- 
nel (p /i)n=i-l-0, the system crosses diabatically the chan- 
nel p+ {fi 0)n=7- Muon transfer is completely negligible 
as the coupling is very small compared with the colli- 
sion energy. The couplings to channels p+ (/iO)„=6 a-nd 
p+{lJ. 0)n=5 are larger as evidenced by avoided crossings. 
The other channels p + {jJ, 0)„<5 are weakly coupled to 
the initial one and they are not expected to be populated 
significantly. 

In Figs. 131 and ^ we present the muon transfer proba- 
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FIG. 4: Multichannel calculated partial and total transfer 
probabilities for the Thomas- Fermi potential. The dashed 
curve is the result of the Eq. 11311 

bilities into the different product channels p+ {iJ- 0)„=5_6 
together with the total transfer probability, for the pure 
coulombic (C) and the Thomas-Fermi (TF) potentials, 
respectively. The two channels are about equally popu- 
lated in both cases and the total transfer probability is 
high (80 %) for energies in the intermediate range (be- 
tween 10"^ and 10^ eV). This can be understood in terms 
of non-adiabatic transitions in a simple 3 channels model 
( {p ^)n=i + O and p -\- {^0)n=6,b)- Using the Landau- 
Zener expressions for the two crossings, a total transfer 
probability of ^ 0.6 between and 10^ eV and a slow 
decrease at higher energies is obtained. This is in quali- 
tative agreement with the fully converged results. 

The main difference between the results for the C and 
TF potentials is at low energies. Whereas for the C po- 
tential the probabilities are constant at low energies, for 
the TF potential there is a typical threshold behavior. 
This is what is expected for short range potentials [Tll| . 
Indeed, in the colinear configuration the C potential be- 
haves asymptotically as 1 /B? (charge-dipole interaction) 
while the TF potential goes to zero much faster. At low 
energies, the amplitude of the wave function for po- 
tentials near the origin is constant, while for short range 
potentials it decreases as energy goes to zero Thus 
the low energy dependence of the transfer probabilities 
is enterily determined by the asymptotic behavior of the 
potentials. This suggests that the dynamics can be de- 
scribed by a two step mechanism: a transmission in the 
incident channel through the long range part of the po- 
tential, followed by non adiabatic transitions in the in- 
teraction region. The charge transfer probability for low 
energies will be then given by 

P{E)^\T{E)\''P^ax (13) 

where T[E) is the transmittance in the incoming channel 
and Pmax the plateau transfer probability. 

The transmittance T{E) can be estimated by a simple 
one channel calculation with an effective potential of the 
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probability. Since the colinear configuration is the most 
favorable configuration, this procedure provides an up- 
per limit of the total transfer probability. The results 
are presented in figure [3 together with those of the co- 
linear calculations. Clearly, the threshold is displaced 
towards higher energies (10~^ eV). Thus we expect that 
for thermal energy collisions (4 10^^ eV) the muon trans- 
fer probability will be of the order of 0.6 or less and will 
increase smoothly with energy in the range 10^'^ — 10^^ 
eV. 



IV. DISCUSSION 



FIG. 5: Total transfer probabilities for the colinear model 
((C): pure Coulombic potential, (TF): Thomas- Fermi model) 
and 3D estimate obtained from Eq. 1131 



form 



^efr(-R) = — " 



3?i^ 
2 m 



(14) 



for R > Ro and V^q{R) = V^q{Ro) for R < Rq. The 
transfer radius has been taken at the avoided crossing 
with the p + (/i 0)„=6 channel. For R < Rq the solution 
of this one-channel problem is given by 



F{R) = T{E) 



sm{KR) 



(15) 



with K = {2mo,p^{E ~ V^Q{Ro))y^^ /n, while for R 
oo we have F{R) = sin(fci? -I- 5)/k^/'^ with k — 
(2mo.ppi?)^/^/?i. Numerical solution of the one chan- 
nel problem thus provides directly T{E). For the pure 
Coulombic potential {Z* — 8) we found T{E) = 1, con- 
sistent with the result of our multichannel calculations 
which gives a constant transfer probability at threshold 
(see Fig. O . This unusual behavior is due to the par- 
ticular 1/R^ dependence of the effective potential [Tl|. 
The shorter range TF potential gives a more standard 
threshold behavior, with a transmittance going to as 
energy decreases. This is apparent in Fig. ^ where the 
total transfer probabilities for the TF potential as a func- 
tion of the energy are presented for both, the multichan- 
nel calculation and the simple estimation using Eq. (|13|l . 
The two calculations agree within 10%. 

This suggests a simple procedure for estimating the 
energy dependence in 3 dimensions (3D). In 3D, the ion- 
dipole interaction perturbes thep/i„=i initial state at sec- 
ond order through a — 1/ i?^ term 8] . So even for the pure 
coulombic case there will be a threshold behavior at low 
energies. We have thus performed a one channel calcu- 
lation of the transmittance T{E) with the 3D functional 
form of the potential. This T{E) was used in Eq. H13|) 
together with Pmax provided by the colinear Coulombic 
potential to yield an approximate 3-dimensional transfer 



In order to make contact with the experiments of 
Werthmiiller ei a/ , we have calculated the muon trans- 
fer rate according to: 



X{E) ^Nva{E) 



(16) 



with N being the number density of liquid hydrogen 
(4.2510^^ cm~^), V = ^2 E/ruo^p^ the relative velocity 
and 



-iE)^y,P{E) 



(17) 



the muon transfer cross section with k — yl^ 



E/h. 



In writing Eq. H17|l . we have assumed only one partial 
wave with zero angular momentum. This is valid for 
collision energies up to 0.1-0.2 eV, corresponding to the 
height of the barrier in the entrance channel for the p 
partial wave. 

In Fig. we present our 3D estimates together with 
the results of the calculations of Sultanov and Adhikari 
based on a two-states approximation of the Fadeev- 
Hahn equations. The agreement is good, our results be- 
ing slightly higher than theirs. This was expected since 
our 3D estimation is based on the colinear calculations 
which provides an upper limit to the transfer probabil- 
ities. In addition, they have consider only the 0/x„=5 
final channel From Eqs. H16|) and H17|l it follows that for 
a constant transfer probability P{E) — 1 the rate should 
decrease as I/VE. This is also plotted in Fig. |S| We 
note that our calculated rates decrease with energy slower 
than those for a constant transfer probability. The rea- 
son is that our calculated probabilities actually increase 
(although slower than 1/y/E) in the energy range con- 
siderde here as can be seen in Fig. |31 

In Fig.Elwe have also represented the rates of the two- 
component model used in Ref. [3| to fit the experimen- 
tally observed X-ray emission from excited muonic oxy- 
gen. While for thermal energies there is an order of mag- 
nitude agreement between our estimated values and the 
two-component model, there is a clear disagreement for 
epithermal energies. The two-component model predicts 
an increase of the rates going from thermal to epithermal 
energies in contradiction with our predictions and those 
of Sultanov and Adhikari (S)]. Moreover, for energies 
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FIG. 6: Estimated transfer rates obtained in this work (3D). 
The dashed curve represents the rate corresponding to a max- 
imum transfer probability P(E) of unity. The points are 
the calculated values of Sultanov and Adhikari _5] using the 
Fadeev formalism. The solid straight lines (2C) are the rates 
determined from fitting the experimental data of Ref. Q with 
the two-component model. 



above 0.1 eV the rate predicted by the two-component 
model is three times larger than the maximum values of 
the rate obtained by assuming a muon transfer probabil- 
ity per coUision of unity (see Fig.O. 



V. CONCLUSIONS 

We have presented colinear calculations of muon trans- 
fer rates between muonic hydrogen and oxygen for rela- 
tive translational energies between 10~^ and 10'^ eV. For 
the lower energies (below 10~^ eV) the de Broglie wave- 
length is much larger than the characteristic distance of 
the potential interaction and the problem corresponds to 
an ultra-cold collision. A simple procedure to estimate 
the energy dependence of the muon transfer rate in three 
dimensions was proposed. 

Our results show that the muon transfer rate decrease 
going from thermal to epithermal energies, in agreement 
with previous theoretical calculations |a] but in contra- 
diction with the assumptions of a two-component model 
used in Ref. |^ to interpret the experimental data on this 
system. 

In the experiment the X-ray emission shows a bi- 
exponential behavior. In the two-component model of 
Werthmiiller et al ;4!], it was assumed that the prompt 
emission is due to muon transfer from epithermal muonic 
hydrogen while the delayed emission, similar to what is 
observed when oxygen is replaced by other atoms, is asso- 
ciated to thermalized muonic hydrogen. Assuming that 
the relative amounts of epithermal and thermal species 
are equal, Werthmiiller et al concluded that the rate of 
muon transfer from epithermal muonic hydrogen to oxy- 
gen should be larger than the one from thermalized ones 



by a factor of almost 4. They suggest that a resonance 
could be responsible for this increase. We have shown 
that for s-waves it is not possible to obtain such an in- 
crease in the transfer rate, even if assuming a transfer 
probability per collision of one. Thus, if resonance effects 
exist they have to be due to the contribution of higher or- 
der partial waves. In an elastic cross section calculation, 
Kravtsov et al pol | have shown that for muonic hydrogen 
colliding with an oxygen nuclei, the contributions of p 
and d waves are small as compared to the s wave in the 
range 0-0.2 eV and that only at ~1.5 eV there is narrow 
resonance involving the d wave. This will thus imply that 
the epithermal muonic hydrogen has much higher trans- 
lational energy than what has been assumed in Ref. Q. 
Anyway, this conclusion is based on a calculation of the 
elastic scattering. It could be of interest to check whether 
it is valid also for the muon transfer process as well. 

In this work the calculations were performed in a re- 
stricted colinear configuration. This approximation is 
known to give poor quantitative results in many cases. 
However for the system considered here, the colinear con- 
figuration is the most favorable for the transfer process. 
Also, since the colinear configuration provides the mini- 
mum energy path for the reactance channel, even for an 
initial non colinear configuration there will be efficient 
orientational effects in particular for low translational en- 
ergies. In addition we have introduce a correction to the 
low energy dependence of the rate in order to take into ac- 
count the correct asymptotic behavior of the potential in 
3-dimensions. Therefore, we believe that although quan- 
titatively the results may change if a full 3-dimensional 
calculation is performed, the conclusions, in particular 
concerning the energy dependence of the muon transfer 
rates, will remain valid. A full 3-dimensional treatment 
is not impossible however, although quite computation- 
nally demanding. It will definitetly be very important 
to further improve our detailed understanding of this re- 
action. It would also be most interesting to study the 
energy dependence of the muon-transfer rates for other 
elements, such as C or S, for which the rates seem to 
behave differently as a function of the energy. This work 
is in progress. 
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